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Quantum entanglement at the ^(3770) and Y(4S) 
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We review results which explicitly depend on the entanglement of neutral meson pairs produced at the -0(3770) 
and T(4S). Time-dependent CP-violation analyses at the B-factories use the flavour-singlet final state at 
the T(4S), but by assuming its quantum-mechanical evolution; Belle on the other hand has tested the time- 
dependent flavour correlation of the B-mesons, comparing predictions of quantum mechanics, spontaneous dis- 
entanglement, and Pompili-Selleri models. At the 0(3770), decay rates are modulated by various combinations 
of the charm mixing parameters: this has been exploited by CLEO-c to provide the first effective constraint 
on the strong-phase difference 5. Finally, the goal of a "model-independent" 03/Dalitz analysis is now within 
reach, using D-mesons from the -0(3770) to constrain the D° K^7r + 7r _ decay amplitude. Manifestly entangled 
events -0(3770) (K^7r + 7r _ )D(K^7r + 7r _ )D, rather than just "CP-tagged" decays, turn out to be crucial. 



1. Introduction 

Einstein- Podolsky- Rosen correlations [1], in the 
form cited by Bohm [2], are one of the most cele- 
brated features of quantum mechanics. For a spin- 
singlet state of photons or particles, 
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measurements on particle 1 (2) are indeterminate, 
but once made they fully determine the result of a 
measurement on particle 2 (1). For the arrangement 
shown in Fig. 1, Bell's theorem [3] (in the CHSH [4] 
form) shows that 

\S\= E(a,b) - E(a,b') +E(a',b) +E(a',b f ) <2 

for any local realistic model; quantum mechanics al- 
lows for \S\ as large as 2\f2 for an optimal choice of 
settings. 1 Results breaching the bound (2) thus rule 
out local realism, even if we subsequently "get be- 
hind" quantum mechanics to a more complete theory: 
the quantum weirdness, or at least this part of it, is a 
phenomenon of nature. 
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Figure 1: Bohm's version of the EPR "experiment". 
From [5]. 



E(a, b) 



-•^Here E is the^ correlation of measurements 

i? ++ (a,6)+i?__ (a,b)-R+- (a,b)-R- + (a,b) 



As shown by Aspect et al. [5], and many times 
since, the bound (2) is broken for photons produced 
in cascade decays, with correlated polarizations: the 
data moreover are consistent with quantum mechan- 
ics. (The Aspect result was S = 2.697 ± 0.015, cf. 
S QM = 2.70 ±0.05.) 

This behaviour is less well-tested for massive sys- 
tems, but is nonetheless routinely used. Decay of a 
vector meson to a pair of neutral pseudoscalars, 
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(|B )|B°)-|B°)|B )), (3) 



R ++ (a,b)+R__(a,b)-R + _(a,b)+R_ + (a,b) 



produces a B-pair entangled in a flavour singlet state: 
the flavours of the individual B-mesons are indetermi- 
nate, but at a given time t, the pair is always B°B°. In 
the flagship measurements of time-dependent CP vio- 
lation at the B-factories, one B-meson B^ AG is recon- 
structed in a state of definite flavor, and the other B GP 
in an eigenstate of CP. The decay rate is then modu- 
lated in the difference of decay times At = t\ — £2, 

-|At|/r B o 

Tcp(At) = — [1 ± {S C p sin (Am At) 

4r B o 

+A CP cos (Am At)}} , (4) 

with one rate for B^ AG (+) and another rate for Br^ AG 
(— ). The coefficients Sep and Aqp in (4) are thus 
CP-violating. 

A complementary measurement, assuming the B- 
physics but testing the time-dependent flavour oscil- 
lation due to the entanglement (3), has been per- 
formed by Belle (Section 2). In the formally equiv- 
alent, but experimentally different setting of ^(3770) 
decays, decay rates of D-meson pairs are modulated 
by the charm mixing parameters (x, y) and the strong 
phase difference 5: CLEO-c has exploited this to de- 
rive a constraint on the latter quantity (Section 3). 
The final state ^(3770) (K° s 7r+7r-) D (K° 7r+7r-) D , 
where the entanglement of the D-mesons is clearly 
manifest, also turns out to be crucial for the "model- 
independent" measurement of the unitarity angle ^3 
using the Dalitz analysis method at the B-factories 
(Section 4). 



Flavor Physics and CP Violation Conference, Taipei, 2008 




Figure 2: An optical analogue of the T(4S) BB EPR 
correlation analysis. 

2. T(45): EPR correlations at Belle 

In the quasi-spin analogy initially expounded for Ja- 
mesons [6], a |B°) corresponds to a spin \^) z particle 
or a photon with vertical polarization (V); a |B °) cor- 
responds to a spin | ty) z particle or a horizontally polar- 
ized photon (H). While optical measurements can be 
made on arbitrary axes a\ ff) + (3\ Jj), only the | ff) and 
| J|) measurements are practical for B-mesons. How- 
ever, over time t the state |B°) evolves to 

X - [{1 + cos(Am d *)}|B°) + {1 - cos(Am d t)}|B °)] , 

making other "measurement axes" accessible. Thus 
for B-pairs produced in the flavour singlet state (3), 
the decay time difference Ara^At plays the role of 
the difference Acj) between polarimeter orientations in 
an optical experiment. An optical analogue of the 
situation at the B-factories is shown in Fig. 2: the 
polarimeters are fixed such that they always perform 
a V/H measurement on the photons, but phase ro- 
tations and (j)2 are inserted between the produc- 
tion point of the entangled photons and their mea- 
surement. 

This setup would seem to allow a full Bell inequality 
test, but in fact the rotations are imposed by 
nature — in the decay times of the particles — rather 
than being subject to the experimenter's choice. Thus 
a true Bell test cannot be performed. 

2.1. The Green Baize Table Conspiracy 
Model 

While there are general arguments for the impos- 
sibility of such a test [7], 2 it is more entertaining 
to consider counter-examples. A breathtakingly un- 
realistic local-realistic model, indistinguishable in its 
predictions from quantum mechanics, was devised for 
this purpose by Bramon, Escribano, and Garbarino [8] 




Figure 3: According to one local- hidden- variable model, 
the Cigarette-Smoking Man takes a close interest in 
T(45) -> BB decays. 



(following Kasday [9]): we expound it here in a form 
inspired by one of the myths of our age. 

Somewhere, there is a wood-panelled room with a 
green baize table, where powerful men meet together, 
smoke, and make conspiracy (Fig. 3). They determine 
world events in detail, including decays T(45) —> BB. 
At the time of each decay, t = 0, four hidden variables 
are set: mesons 1 and 2 are each given a piece of paper 
(as it were) bearing the quantities (ti, /i) and (£2, /b)- 
These act locally: meson i decays at a time t = into 
final state / = according to the values written on 
its piece of paper. Now as part of the conspiracy, in 
order to deceive the world, (£1, /1, £2? /b) are chosen 
randomly according to quantum mechanical rules. 

The phenomena in this model are indistinguishable 
from those in quantum mechanics, even though the 
individual particles have a definite state at all times: 
the model is local, and (if only in this sense) realis- 
tic. 3 Another explicit counter-example has since been 
constructed by Santos [11]. 



2.2. QM versus specific LR models 

The best that can be done is thus to test the 
data in detail against quantum mechanics, and the 
predictions of specific local realistic models. Rea- 
sonable, non-conspiratorial models can then poten- 
tially be ruled out. The experimental quantity of 
choice is the time-dependent asymmetry of decay rates 
to opposite-flavour (OF) and same-flavor (SF) final 



2 It turns out that even in the case of active flavor mea- 
surement (rather than letting the mesons decay), a Bell in- 
equality test cannot be performed, because the decay of the 
B-mesons is too rapid compared to their flavour oscillation, i.e. 
Xd = Arrid/^d — 0.77 is too small. 



3 Optical experiments do not share this vulnerability. For ex- 
ample, in the modification of the Aspect experiment reported by 
Weihs et al. [10], the effective polarimeter orientation is changed 
while the photons are in flight, according to randomly gener- 
ated numbers. In this case, no conspiracy can fix the results to 
conform to a particular pattern in A(p = Ani^At. 
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states: 

A(tl ' t2, -iJ0F(tl,t2)+i25F(tl,*2) (6J 

= cos(Am d At) (7) 

for quantum mechanics. The dependence on At = 
ti — ^2 alone, a consequence of the entanglement of 
the state, is distinctive. It's therefore useful to con- 
sider asymmetries as functions not of (ti,^), but of 
At and t m i n = min(ti,t2). Several possibilities for 
A(At, t m in) are shown in Fig. 4. As a limiting case, 
under spontaneous disentanglement (SD) to a B°B 
pair (with definite flavour) immediately after T(4S) 
decay, the two mesons undergo independent flavor os- 
cillations, with 

Asd = cos(Ara^i) cos(Ara d t 2 ), (8) 

taking the complicated form in the figure when plot- 
ted on (At, t m i n ). More seriously, one can consider 
the class of models obeying the assumptions of Pom- 
pili and Selleri [12]: QM-like states for the individual 
mesons with stable mass; 100% flavour correlations, 
to reproduce the QM behaviour as closely as possible; 
and the constraint that QM predictions for uncorre- 
cted B-mesons should also be preserved. Asymme- 
tries must then lie between the bounds 

Apf 1 = 1 - min(2 + *, 2 - *), where 
\I/ = {1 + cos(Am^At)} cos(Am^t m i n ) 

- sin(Ara d At) sin(Am d t min ) (9) 

and 

Ap| x = 1 - | {1 - cos(Am d At)} cos(Am d t min ) 

+ sin( Am d At) sin( Am d t min ) | (10) 

shown in the figure. 

In the current state of the art, the discrimination 
power shown in Fig. 4 is not fully realised, as only At 
is measured: the expressions (8)-(10) must in effect 
be integrated over t m i n . It turns out that it is still 
possible to exclude the local-realistic models shown. 

2.3. The 2007 Belle result 

The Belle analysis [13] is based on a 152 x 10 6 BB 
data sample, and uses techniques established for time- 
dependent CP-violation analyses: one B-meson is re- 
constructed in a flavour-tagging mode B° — > D*~^ + ^, 
while the other B-flavour is tagged using a lepton; a 
consistency check with other flavour- tagging informa- 
tion is imposed to maintain purity. A sample of 8565 
such events is found: 6718 opposite-flavour, and 1847 
same-flavour pairs, divided into 11 bins in At. 




ASYMMETRY A(At,t min ) for POMPILI-SELLERI (min) 




Figure 4: For ideal measurement at T(4S) — > BB: de- 
cay rate asymmetry A — (Rof — Rsf) /(Rof + Rsf) 
as a function of (At, t m i n ), (top) for quantum mechan- 
ics; (second) for spontaneous disentanglement; and for the 
Pompili-Selleri class of models, showing (third) minimum, 
and (bottom) maximum values. 
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At [ps] At [ps] At [ps] 

Figure 5: From [13]: Asymmetries as a function of At in Belle data, compared to (left) quantum mechanics, (middle) 
spontaneous disentanglement, and (right) Pompili-Selleri models. The lower plots show time-dependent flavour asym- 
metry (crosses) and the results of weighted least-squares fits to each model (rectangles, showing ±1<t errors on Awj). 
Upper plots show differences A = Adata — ^modei in each bin, divided by the total experimental error <j to t- Bins where 
Apg n < Adata < Apg x have been assigned a null deviation. 



Backgrounds are subtracted from the OF and SF 
samples separately, in At bins: fake D*, using side- 
bands (126 ± 6 OF, 54 ± 4 SF); bad D* - i combina- 
tions, also from data (78 ± 9 OF, 236 ± 15 SF); and 

B + — > D iv events, estimated from Monte Carlo 
and the cosb,d*£ distribution in data (254 OF and 1.5 
SF [both ±6%]). This process is crucial, as the back- 
grounds produce a time-structured difference in the 
asymmetry. After a (1.5 ±0.5)% correction for mistag- 
ging, deconvolution is performed (using DSVD, on 
the SF and OF samples separately) to remove decay- 
vertex-resolution, efficiency, and other remaining ef- 
fects. (Potential bias against any of the models is 
explicitly studied and subtracted, with a systematic 
error assigned.) Finally, a simple exponential fit to 
the resulting histogram is performed, to extract the B- 
lifetime as a check: r^o = (1.532 ± 0.017) ps is found, 
cf. the (1.530 ± 0.009) ps world average [14]. 

A comparison between the data, and fits to the 
predictions of each model, is shown in Fig. 5. In 
each case, the parameter Am^ is floated, subject 
to the constraint of the world average value, with 
Belle and BaBar contributions excluded: (0.496 ± 
0.014) ps -1 [15]. 4 Quantum mechanical predictions 
fit the data well (x 2 / n dof — 5/11), while spontaneous 
disentaglement is disfavoured at 13a (x 2 = 174). A 
fit to (1 - C b ob o)A QM (At) + C b ob o Ast>(A£) finds an 
"SD fraction" of C b obo = 0.029 ± 0.057, consistent 
with zero. 5 The entire class of local realistic models 
satisfying the minimal Pompili-Selleri assumptions is 
disfavoured at 5.1<j (x 2 =31). 



4 The Belle and BaBar results are from fits assuming QM 
time evolution, and so cannot be used. They otherwise domi- 
nate the average, improving its precision by a factor of three. 

5 In a decoherence model, this is equivalent to multiplying the 
interference term in the B°-B basis by a factor (1— CbOb°) [1^1- 



3. ^(3770): (x, y, S) and rates at CLEO-c 

Formally, the situation at the ^(3770) is the same 
as that atthe T(45): the decay e+e~ -> ^(3770) -> 
(|D°)|D°) - |D°)|D )) leads to an entangled final 

state equivalent to (3). However there are differences 
in practice: mixing is a percent-level effect in D-decay 
amplitudes, and CP violation is suppressed orders of 
magnitude further. So while the principal use of the 
state (3) at the T(4*S) is study of CP violation, an 
obvious use of D-meson pairs produced at the ^(3770) 
is CP tagging. For example, decays to two CP-even 
(or two CP-odd) eigenstates don't occur. 

If we consider decays ^(3770) — > 
(K - 7t + )d(K - 7t + )d, however, the situation is not 
so straightforward. Relative to production from a 
pair of D° mesons, the rate is suppressed by the 
mixing rate Rm = \( x<2 + y 2 )'-, by contrast, the rate 
for uncorrelated DD decays to this final state is 
suppressed by only the "wrong-sign" rate Rws, it 
is forty times larger. There are thus nontrivial effects 
due to the coherence of the state produced at the 
^(3770). Currently the most systematic treatment 
is [17], following earlier work by [18] and others. 

3.1. D° K° S L 7T 

A simple example is the study of D° — > K^7r° by 
CLEO-c [19], where the D° K°tt° decay is recov- 
ered using the distribution of missing-mass-squared 
M^ iss in events tagged by a fully-reconstructed D 
decay. In prat ice there are three distinct samples, cor- 
responding to the three tagging modes D° ^ K + 7r~, 
K+tt-tt , and K+7r-7r-7r+. A D° K°tt° branch- 
ing fraction calculation using tagging mode / in fact 
determines 

/ 2r f cos5 f +y \ 
Bk °^° V + l + Rwsj ) ' 
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Table I From [20]. Effective branching fractions (upper 
section) for D° decay modes i, divided by and (lower 
section) for D°D decay to modes divided by BiBj, 

in V(3770) -> ^ (|D°)|D°) - |D°)|D )) data. S + (5_) 
denotes a CP-even (CP-odd) eigenstate, and e~ a semilep- 
tonic final state X + e~u e - Quantities are shown to leading 
order in the mixing parameters (x,y) and the wrong-sign 
rate Rws- Rm — \{x 2 + y 2 ) is the mixing rate. 



Mode 


Correlated 


Uncorrelated 


K~7T+ 








2 


2 


S- 


2 


2 




Rm 


i?WS 




[(1 + Rws) 2 


1 + Rws 




— 4r cos 5 (r cos 6 + y)] 






1 + Rws + 2r cos 8 + y 


l + Rws 




1 + Rws — 2r cos S — y 


1 + Rws 




1 — ry cos 6 — rx sin 5 


1 







1 


S-,S- 





1 




4 


2 




1 + 3/ 


1 


S-,e~ 


1-2/ 


1 



3950108-009 




1 2 -0.010 0.000 0.010 0.020 
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Figure 6: From [20]: Constraints on the strong phase 
difference 5, after combining results from CLEO-c's en- 
tangled source of D-mesons with external branching frac- 
tion and mixing parameter measurements. Likelihood 
curves (including statistical and systematic uncertainties) 
are shown for (a) cos£, (b) xsin£, and (c) 5; (d) shows 
contours in units of yj Ax 2 on (cos 5, x sin 6). 



where the amplitude ratio r/ and strong phase differ- 
ence S f (r f e~ i6 f = (/|D°)/(/|D )) and the wrong- 
sign rate Rwsj are mode-dependent. 

The method is to use the equivalent measurement 
for D° — > K^7T°, an an untagged measurement of the 
branching fraction for that mode, to determine the 
product Cf = (2rf cos5f+y)/(l + R W sj). (A tagged 
analysis in this mode finds an effective branching frac- 
tion S k Otto(1 - C f ).) The true D° -> K°7r° branching 
fraction can then be measured for each tagged sample: 
averaging over them, CLEO-c find 

= (0.998 ± 0.049 ± 0.030 ± 0.038)%, (11) 

and an asymmetry 

By-0 — o 

= 0.108 ± 0.025 ± 0.024, (12) 

consistent with the value 2tan#c = 0.109 ± 0.001 
expected if symmetry under the U-spin subgroup of 
SU(3) is imposed. 

3.2. Charm mixing and S 

In the general case, effective branching fractions 
depend on the mixing parameters (x,y), and (for 
the K~7T + final state) the strong phase difference 
5 ee 5k<k- Correction factors are summarized in Ta- 
ble I. Whereas in the D° — > K^7r° analysis the correc- 
tions were a complication that needed to be taken into 



account, the CLEO-c analysis reported in [20] uses a 
suite of measurements, and their varying dependence 
on (x, y) and 5, to constrain the mixing parameters, 
and in particular S. 

In the first case a least-squares fit is performed 
to the yields in eight hadronic final states (K~7r + , 
K+7T-, K+K-, 7T+7T-, K°7r°7r°, K°tt°, K°t?, and 
K° s u), and 43 "double-tagged" final states (24 fully- 
reconstructed, 14 including a semileptonic decay, and 
5 including K^7r°), together with external results on 
seven branching fractions (CP eigenstates, and K~7r + , 
with correlations taken into account). The result finds 
the lifetime difference parameter y with large uncer- 
tainty, and as a result the combination #sin<$ is un- 
constrained. 

An extended fit, including measurements of mixing- 
related quantities (y, x, r 2 , y' , (x f ) 2 ) by other exper- 
iments, is therefore performed: results are shown in 
Fig. 6. The fit finds 

cos 5 = 1.10 + 0.35 + 0.07 

x sin S = (4Atli ± 2 - 9 ) x 10 " 3 > ( 13 ) 

and after minimising on the physical surface 
(cos 5, sin (5), 

£=(221^)°, with 

S e [-7°, +61°] at 95% confidence, (14) 

the first effective constraint on this quantity. To- 
gether with the external mixing measurements (with- 
out which the analysis is not possible), the precision 
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on cos 5 is driven by the yields for the eight hadronic 
final states, and double-tag yields {K-zr, S±} including 
CP-eigenstates S±. These results are for a 281 pb -1 
sample: better precision will be possible with the final 
CLEO-c dataset. 



4. 3 /Dalitz: ^(3770) rescues the T(45) 



Dalitz analyses of B d 



DK ± and D*K ± , with 
D — > K^7r + 7r _ and K^K + K~, are currently the most 
sensitive probe of the unitarity angle 03 (also known, 
in the least interesting of the disagreements between 
the B-factories, as 7). The state of play in these 
important analyses was shown in Anton Poluektov's 
talk on Monday [21]; Belle's preliminary update [22], 
shown in Fig. 7, finds </> 3 = (76±J§ ± 4 ± 9)°; BaBar's 
new publication [23], shown in Fig. 8, finds 03 = 
(76 ± 22 ± 5 ± 5)°. In both cases, the model error 
(shown last) is already uncomfortably large. 

The future of this measurement is the so-called 
model-independent approach [24] ; the feasibility study 
for the method has recently been updated [25]. Rather 
than relying directly on (say) an isobar model to deter- 
mine the D — > K^7r + 7r - decay amplitude, the phase 
difference parameters 



c = cos((5£>(m^_, m 2 _) — ^(m 2 ,, m\)) 

s = sin(fe(m^_, m 2 _) — fe(m?_, m+)) (15) 

are extracted from ^(3770) data, exploiting the cor- 
relations in the final state to measure the Dalitz plot, 
not for D° or D°, but for CP-tagged D-mesons. The 
challenge is to cope effectively with the limited data 
samples available (or foreseen). 

An advance has been the nontrivial binning shown 
in Fig. 9, which is unifotm in Afe| moder However, 
results prove to be biased for finite Dcp — > K^7r + 7r~ 
sample sizes: tests where events are generated in one 
model, and reconstructed using another, find a return 
of model dependence, as shown in Fig. 10. It arises 
because for each bin, of the parameters in (15), only 
Ci is reconstructed; Si is recovered by a cf + s 2 — 1 
constraint. A change of model results in a shift of the 
Sd region sampled by a given bin, introducing a bias 
in Si via the constraint. 

The study [25] finds that if {q,s^} are determined 
from ^(3770) — > (K^7T7r)D(K^7r7r)D events, the out- 
come is unbiassed for finite data samples: a change 
of model can degrade the sensitivity (Fig. 10), but 
not introduce a bias. The importance of the entan- 
gled final state at the ^(3770) thus goes beyond "CP- 
tagging" : the additional correlations in the final state 
(K^7T7r)D(K^7r7r)D are crucial. 
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Figure 7: From [22]: (preliminary) updated Belle results. 
Dalitz distributions for D — > K S 7T + 7r~ decays from (top) 
B ± -h. DK ± and (bottom) B ± -> D*K ± , for (left) B + and 
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Figure 8: From [23]: a selection of (published) BaBar re- 
sults. Dalitz distributions of (Upper) D — > K£7r + 7r _ 
for (a) B" DK" and (b) B + -> DK+; and (lower) 
D° -> K^K+KT for (a) B" -> DK~ and (b) B+ -> DK+. 




m 2 . (GeV7c*) 



Figure 9: From [25]: Uniform binning in ASD\ model on 
the D° -> K^tt+tt" Dalitz plot; see the text. 
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Figure 10: From [25], showing results of toy Monte Carlo 
simulation of the 03/Dalitz analysis. (Top:) using Dcp 
data; the histogram shows fit results where the same D° 
decay amplitude is used for generation and binning, and 
the points with error bars show the case with different 
amplitudes. (Middle:) as above, but using -0(3770) — > 
(K° S 7r + 7T~)r>(K° S 7T + 7T~)r> data. (Bottom:) coefficients 
(ci,Si) obtained from the eight Afo| modeZ bins of Fig. 9, 
where (left) the same amplitude, and (right) different am- 
plitudes are used for event generation and binning. 



5. Summary 

Analyses which explicitly depend on the entangle- 
ment of neutral meson pairs are becoming important 
in this field. The effect allows the model-dependence 
in 03/Dalitz analyses to be lifted: CLEO-c data will 
already be necessary to enable full use to be made of 
final B-factory results; BESIII results will be needed 
to exploit the data from a sup er-B /flavor factory. En- 
tanglement in ^(3770) — > D°D modulates tagged de- 
cay rates in a way that must be taken into account for 
D° — » K^7T° measurement, and that has now enabled 
the first effective constraint, 5ktt = (22^2^11)°, on 
the strong phase difference in D° — > K + 7r~. And at 
the T(45), entanglement, used routinely in B-factory 
measurements, has now been tested, even though a 
Bell inequality analysis cannot be performed. A con- 
straint on the decoherent fraction Cb°b° = 0-029 ± 



0.057 is found; as for "realistic" local-realistic models, 
the class obeying the Pompili-Selleri assumptions has 
been ruled out at 5.1<j, the first such constraint from 
data. 
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